Free vortex and vortex-pair lifetimes in classical two-dimensional easy-plane magnets 
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We report numerical simulation results for free-vortex lifetimes in the critical region of classical 
two-dimensional easy-plane ferro- and antiferromagnets having three-component order parameters. 
The fluctuations in the vortex number density in a spin dynamics simulation were used to estimate 
the lifetimes. The observed lifetimes are of the same order of magnitude but smaller than the 
characteristic time-scale above which a phenomenological ideal vortex-gas theory that has been used 
to account for the central peak in the dynamic structure factor S"" (q, uj) is expected to be valid. 
For strong anisotropy, where the vortices are in-plane, the free vortex lifetimes for ferromagnets 
and antiferromagnets are the same, while for weak anisotropy, where the vortices have nonzero 
out-of-easy-plane components, the lifetimes in antiferromagnets are smaller than in ferromagnets. 
The dependence of the free-vortex and total vortex densities on the size dependent correlation 
length in the critical region is examined. We also determined the lifetimes of vortex-antivortex 
pairs for T = Tkt a-nd well below Tkt- The observed time-scales are very short, and the observed 
pair densities are extremely small. These results suggest that pair creation and annihilation are 
unlikely to play any role in the central peak in S^^{(\,ui) observed in computer simulations for the 
ferromagnetic model for T < Tkt- 

PACS numbers; 75.10.Hk, 75.30. Ds, 75.40. Gb, 75.40.Mg 



I. INTRODUCTION 

The interpretation of the central peak in the dy- 
namic structure factor ^""(q, w) in easy-plane layered 
ferromagnets (FM) and antiferromagnets (AFM) is cur- 
rently based on the jphenomenological theory developed 
by Mertens et aluU The central peak is observed for 
T > Tkt (Tkt is the Kosterlitz-Thouless transition tem- 
perature) and the theory accounts for it within the frame 
of a dilute gas of free vortices effectively assuming infi- 
nite free- vortex lifetime. It is of importance to determine 
the free-vortex lifetime to find out the interval of appli- 
cability of the theory as well as to understand the time 
scale of the processes which account for the finite life- 
time. We have developed a method to calculate the vor- 
tex lifetime and we present here our results obtained from 
combined cluster Monte Carlo (MC) and spin dynamics 
simulations. We believe that these are the first reported 
free- vortex lifetimes in AFM for the system Hamiltonian 
we have considered. 

The intensity and width of the central peak predicted 
by the dilute- vortex-gas theory depend on the free vortex 
density Ny and the rms vortex velocity u. A free vortex is 
assumed to exist longer than a characteristic observation 
time 
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nu, 



(1) 



where ^ is the correlation length. It is assumed that the 
vortices have a Maxwellian velocity distribution. If this 
theory is valid, it could, in principle, allow one to deter- 
mine the average free vortex velocity u and the correla- 



tion length ^ at a given temperature by experimentally 
measuring the shape of the central peak, knowing the 
strength of the Heisenberg exchange interaction for the 
magnetic material. 

The FM and AFM static properties on a square lat- 
tice (or other bi-partite lattice) are identical (the spins 
on one of the sublattices in the AFM are inverted when 
compared with the spins in the FM) while their dynamic 
properties differ.S In particular, the out-of-plane tilting 
of spins near the core of a vortes— produces a nonzero 
topological charge (or gyrovectorBO) for vortices in the 
FM model. The gyrovector plays an important role in 
determining how the motion of an individual vortex is 
influenced by effective fields due to neighboring vortices, 
and it appears together with an effective mass in a col- 
lective coordinate equation that describes the motion of 
a vortex center O However, the gyrovector is always zero 
for vortices in the AFM model, and, the effective mass 
for vortices in the AFM naadel is much smaller than for 
vortices in the FM model.E3 For these reasons we expect 
that the FM and AFM models could have quite different 
vortex lifetimes, and therefore contrast the lifetimes mea- 
sured for the AFM model wittuthose we have determined 
previously for the FM modeLEj 



II. THE MODEL 

The Hamiltonian of the model is 

H = -JJ2{S^S^ + SfS^ + XS^S]), (2) 
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where the classical spins Si = {Sf , Sf , S'f ) are vectors on 
the unit sphere S"^, the sum is over nearest-neighbor sites 
of a square lattice, the easy-plane anisotropy parameter 
A varies in the interval < A < 1, and J > for FM and 
J < for AFM. We assume that J has energy units, time 
is measured in units of h/ J and temperature in J/ks] the 
spins are dimensionless. 

The static critical behavior of this model is well- 
described by the extensively studied classical XY model 
(the two-dimensional 0{2) model) but because of the pos- 
sibility for spin fluctuations in one more degree of free- 
dom (the out-of-plane fluctuations) the properties of 
the XY model will be observed at lower temperatures 
here. 

Long-range order campjt exist for any .Mnzero tem- 
perature in these modelsjld but Berezinskiilla showed the 
existence of a phase transition which was understood 
by Kosterlitz and Thoulessll3 to occur through unbind- 
ing of pairs of nonlinear topological excijtations-vortices. 
The critical properties of the XY modeO as well as the 
low-temperature phase are of considerable interest to the 
physics of low-dimensional magnetism because there are 
a number of materials which are effectively considered 
as two-dimensional magnetic systemsta and the Hamil- 
tonian (||) is the simplest model to study their proper- 
ties. These are layered crystals with intralayer exchange 
interaction much greater than their interlayer exchange 
(more than two or three orders of maepittijde greater) and 
there are observed both ferromagnetst3"E3 (e.g. K2CUF4, 
(CH3NH3)2CuCU rW2H5NH3)2CuCl4, Rb2CrCl4) and 
antiferromagnet£3 E3 

(like BaNi2(P04)2, BaCo2(As04)2). The intralayer ex- 
change interaction is more than 10^'^ orders of magni- 
tude smaller than the intralayer exchange in the graphite- 
intercalated CQCI2 which has been studied in detail fa^ 
Wiesler et alS3 There are magnetic lipid monolayers ,c3 
e.g. the compound Mn(Ci8H3502)2, which are true two- 
dimensional magnets. 

When the anisotropy parameter A is Less than a criti- 
cal value, Ac which is lattice dependent,c3 only in-plane 
static vortex spin configurations ((5*^) = 0) are sta- 
ble, while for A > Ac only static vortices with nonzero 
out-of-plane ((5^) ^ 0) spin components are stable. 
The criticaJ_anisotropy parameter Ac has recently been 
determinedtj with good precision and for the square lat- 
tice Ac « 0.7034. 



III. DYNAMICS 

The differences in dynamical properties of FM vs. 
AFM not only appear in the vortex properties, but also in 
the spinwave-vortex interaction and the spinwaves alone. 
The spin-wave excitations have a single branch for FM 
while there are two branches for AFM due to the two 
different spin sublattices. The spin waves are either in- 
plane or out-of-plane depending on the value of A. Ivanov 



et aim found a localized mode for the out-of-plane vor- 
tex in AFM which appears in the, gap between the two 
spin-wave branches. Wysin et alx3 have also shown that 
even for the in-plane vortices in the AFM model, this 
localized mode is still present. For the FM model, in 
contrast, only quasi- local spinwave modes appear on a 
vortex (where the component of the mode is localized 
near the vortex core while the , components of the 
mode are extended). _ _ 

In the ideal vortex-gas theory,ErLl the motion of free 
vortices in a temperature interval just above Tkt is as- 
sumed ballistic and leads to central peaks in both in- 
plane S^^{q,uj) [S^^ — due to the symmetry of the 
Hamiltonian) and out-of-plane S^^{q,Lu) dynamic struc- 
ture factors, where S"'°'{q,u;) (a = x,y,z) is the space- 
time Fourier transformation of the correlation function 
(S'"(r,t)5"(O,0)). The central peak (CP) of the in-plane 
dynamic structure factor is predicted to have a squared 
Lorentzian shape 



27T [c^2+^2(l + (^^*)2)]2' 



(3) 



where 7 is defined in (|i|). The CP is located at q* = (0, 0) 
for FM (q* = q) and at the Bragg point K = (tt, tt) for 
AFM (q* = K - q). A weaker CP is predicted for the 
out-of-plane dynamic structure factor for any A, due to 
correlations caused by vortex motion. It is located at q = 
(0, 0) for both FM and AFM. Its intensity is proportional 
to the free vortex density Ny and has Gaussian shape 



exp 
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(4) 



where fx = 1 - A for FM and 1 + A for AFM. For 
A > Ac, the out-of-plane vortex spin asymptotic behav- 
ior is known S'^(r) ~ -^/r^/r exp(— r/ri,), where = 

•\/A/(l — A) /2 is considered as the vortex core radius and 
r is the distance from the vortex center. This asympt. 
totic forrojis used to calculate S^^{q,uj) for both FMEI 
and AFMQ'lI leading to an additional CP only for A > Ac 
with Gaussian shape 



5'^^((7* 



NyU 



(5) 



and again q* = q for FM and q* = K — q for AFM. 
By measuring the width and the integrated intensity of 
the CP one can determine u and ^ and compare with 
independent theoretical results. We compare the charac- 
teristic times 1/7 obtained in this way with our results 
on the vortex lifetime. 

Thus, we expect that vortices in FM and AFM for A < 
Ac have quite similar dynamics (and similar lifetimes), 
which is reflected particularly in their similar in-plane 
correlations. However, the nonzero gyrovector for FM 
vortices is present for A > Ac. This is an additional 
term in the vortex equation of motion for the FM case, 
which is absent for AFMs. Furthermore, the effective 
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vortex mass of AFM vortices-jhas been estimated to be 
smaller than for FM vorticeal3 only for A > Ac- Because 
of these differences, we calculate the vortex lifetime for 
both A < Ac and for A > Ac, with the expectation that 
larger differences should occur in the, latter case. 

In recent numerical simulations,tJ CP has been ob- 
served in S''^{q,Lu) for T < Tkt and A = 0. It is not 
clear currently what causes the CP for T < Tkt- For 
these temperatures, the dominant excitations are spin 
waves and the vortices are bound in pairs with oppo- 
site vorticities. The vortex pairs should only renormalize 
the shape of the spin- wave peaks and should not lead to 
additional peaks as do the free vortices above Tkt- Nev- 
ertheless, we have calculated and checked the lifetime of 
vortex pairs below Tkt- We found that the observed pair 
lifetime is a very short time scale to be the relevant effect 
for the observed CP at these temperatures. 



The move is microcanonical and it also reduces criti- 
cal slowing down even though it is a local move. The 
Metropolis algorithm is needed to satisfy ergodicity of 
the MCS for the Hamiltonian we study. The first 1000 
to 5000 MCS were used for equilibration of the system 
and IC were written after each bin consisting of 2000 to 
10000 MCS. Between 25 to 100 IC were produced for the 
spin dynamics simulations. The calculation of the cor- 
relation length and vortex density is from averages over 
20 or 25 bins with measurements after each MCS and 
all estimated statistical errors are equal to the standard 
deviation of the bin averages. 

The corpeJation length is determined by 
measuringO^m 



a = -v/G(0)/G(q)-l, (9) 



IV. THE SIMULATION 

We simulate a system of classical spins (three- 
component unit vectors) on a square lattice with periodic 
boundary conditions and L x L number of sites, for sizes 
L < 256. First, we run Monte Carlo (MC) simulations to 
obtain initial equilibrium configurations (IC) and to cal- 
culate the correlation length and vortex densities. The 
obtained IC's are then evolved in time using a fourth or- 
der Runge-Kutta methocB to solve the Landau-Lifshitz 
equations of motion 

where 

U(i)) 

X, y, and z are unit vectors along the coordinate system 
axes. The sum is over the nearest neighbor sites of the 
site i. 

We used a combination of Wolff's one cluster 
algorithrrO (for updating of the in-plape spin compo- 
nents Sf and Sf) -Metropolis' algorithmJ^J and the over- 
relaxed algorithmE3 to update the spins in the MC part of 
the simulations. Each Monte Carlo step (MCS) through 
the system consists of three single-cluster updates, three 
Metropolis sweeps, and three over-relaxed sweeps. The 
use of Wolff's algorithm is essential for reducing criti- 
cal slowing down. The over-relaxed move on a classical 
spin consists of rotating the spin around the direction of 
the local field due to its nearest neighbor spins (assum- 
ing no external magnetic field) at an angle of tt, i.e. if 
h = hioc/\hioc\, then 

Snew = 2{h ■ Sold)i^~ Sold- (8) 



G(q) = 




where ~ (5*^^, Sf ) is the in-plane spin vector. The 
wave- vector q = (27r/L, 0) is used, and is the radius- 
vector of the ith lattice site. The correlation length 
£,L obtained from Eq. (||) approaches the exact value as 
0{q*) approaches, zero when increasing the linear size 
of the system L£3'c3 The £aite-size effects for the clas- 
sical XY model are small,liil as a rule of thumb, when 
T/a > 6. 

Anqtljicp, method used to determine the correlation 
lengthliil"E3 in the XY model from simulations on a lat- 
tice is from a fit of the zero spatial momentum two-point 
correlation function to 

T{y,) K^.-y.) ■ Sf.„o) ) - Ccosh((L/2 - y,) /O , 

(11) 

the summation is over the x-coordinates of the lattice 
sites. The advantage of using Eq. (^) to calculate the 
correlation length is that it does not require fitting. We 
have also calculated 7{yi) for several temperature val- 
ues and the comparison betweep-|the correlation lengths 
obtained from Eq. (j|) and froma 

C(2/,) = in(r(y.)/r(y.-i)), (12) 

at saturation (neglecting the contributions due to the 
periodic boundary conditions in Eq. (|T^)), shows that 
the two results agree within 2.5%. 

V. RESULTS 
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A. Free Vortices and Correlation Length 

We consider in this section how free vortices are de- 
termined and the associated difficulties when apphed to 
simulations on a lattice. For the calculation of the vof-i 
tex lifetime we use a method to locate "free" vorticesEJ 
which is not directly related to the correlation length and 
therefore it is very important to understand how it relates 
to a method directly using ^(T) as the length scale to 
distinguish free from bound vortices. This requires the 
calculation of the correlation length for different values 
of T and A which we present below. The results for the 
correlation length by themselves should be of general in- 
terest, in particular, the values of ^(T) for A = 0.9 may 
not be reported elsewhere. ■— , 

According to the Kosterlitz and Thouless theoryjlj at 
Tkt vortex-antivortex pairs start to unbind and free vor- 
tices are formed for T > Tkt- Treating the core ra- 
dius of a vortex as a variational parameter and minimiz- 
ing the vortex energy with respect to it, one can show 
that the vortex core radius is proportional to the corre- 
lation length ^(T). Thus, ^(T) sets a length scale be- 
low which vortex-antivortex pairs can still be considered 
bound while above it the vortices are free. The average 
free vortex density becomes 

n4T) ^ lie{T), (13) 

and there are arguments^ that the exact dependence is 
n.(T) = l/{^e{T)). 

This suggests the following approach in a search for free 
vortices on a lattice. A vortex is considered free, provided 
the minimum distance from its center to a nearest neigh- 
bor vortex is greater than ^(T). If the distance is less 
than £,{T), the corresponding two vortices are marked as 
bound. 

There are two related problems which render this ap- 
proach too CPU intensive currently. First, the simu- 
lations are on a lattice and a vortex center can be at 
any point inside a plaquette (1x1 square on the lattice 
with spins located at its corners), where thej-uprtex has 
been found. We use Tobochnik and ChesterEj method, 
which accounts correctly for measuring spin in-plane an- 
gles mod(27r), to determine which plaqucttcs contain vor- 
tices. Fitting to the known vortex solution for the spin 
Si = {cos 9i cos (pi, cos 9i sin (j)i, sin^^), with in-plane an- 
gle 

0,=qtan-if^i^), (14) 

where g = ±1 is the vorticity and (x^, y^) are the co- 
ordinates of the center of the vortex (treated as fitting 
parameters), allows determination of the vortex center 
location but creates the second problem. The vortex 
positions and their free/bound status have to be deter- 
mined at each time step during the spin dynamics part 
of the simulations in order to determine the free vortex 



lifetime. The fitting procedure for precise vortex posi- 
tioning increases the CPU time by more than one order 
of magnitude (and close to two) which makes the simula- 
tions impractically long for large enough lattice sizes (in 
order to have negligible finite size effects). In addition, 
using the correlation length to determine the free vortices 
for S,{T) > 5.0 gives approximately no free vortices. For 
temperatures such that S,{T) > 5.0 the vortices are still 
distributed in clusters and the correlation length is usu- 
ally greater than the minimum distance between nearest 
neighbors of vortices. |— . 

Therefore, we followed our previuos approachcJ for lo- 
cating free vortices. We take the vortex positions as the 
centers of the plaquettes, which eliminates the most CPU 
intensive task of precisely fitting their positions. We also 
use a fixed length scale to determine the free/bound vor- 
tex status in the temperature interval we study. This 
length scale is set equal to the next nearest neighbor dis- 
tance, i.e. \/2 lattice constants on the square lattice. A 
bound vortex will have a nearest neighbor vortex at a 
distance of one lattice constant or a/2 while a free vortex 
may have its nearest neighbor vortex at a distance greater 
or equal to 2 lattice constants. This method to deter- 
mine the free vortices has already been used by Gupta 
and BailliecJ to measure the vortex density. 

We determined the correlation lengths and vortex den- 
sities for the temperature intervals we studied and two 
values of the easy-plane anisotropy parameter (A — 0.0 
and A = 0.9) in order to estimate the applicability of 
the fixed-length (\/2) approach to determine the free 
vortices. With the critical temperatures estimated as 
rKT(A = 0.0) « 0.70 and rKT(A = 0.9) « 0.63, we 
studied the range 0.75 < T < 1.1 for A = 0.0 and 
0.7 < T < 1.05 for A = 0.9 at intervals of AT = 0.05. 
The critical temperatuce was estimated using the reduced 
forth-order cumulant.Cj For T lower than these intervals, 
the free-vortex density is too low to obtain good statis- 
tics. The results for the correlation length are presented 
in Table | for A = 0.0 and A = 0.9. We expect that 
these values of have reached saturation and should be 
approximately equal to the values of ^ in the thermody- 
namic limit with the possible exception of the data point 
at T = 0.75 for A = 0.0 and at T = 0.7 for lambda = 0.9. 

These values of the correlation length suggest that the 
use of fixed length equal to \/2 to determine the free 
vortices, instead of the correlation length, will overesti- 
mate their number for 0.75 < T < 0.85, A = 0.0 and 
0.7 < T < 0.8, A = 0.9. In these temperature intervals, 
a vortex has its nearest neighbor vortex at a distance 
smaller than ^(T), for most of the vortices, and the re- 
sulting free vortex density will be one or more orders of 
magnitude less than the theoretical result Eq. (p^. For 
0.9 < T < 1.1, A = 0.0 and 0.85 <T < 1.05, A = 0.9, the 
average number of free vortices should be approximately 
the same when measured by using the correlation length 
or the \/2 fixed-length scale. The overestimation of the 
free vortex number very close to Tkt will also lead to 
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overestimation of the free vortex hfetime for these values 
of T (see the next subsection). However, this does not 
change the conclusions regarding the magnitude of the 
lifetime with respect to the ideal vortex gas theory. 

We also calculated the dependence of the average free 
and total vortex densities on 1 to study their spatial 
distribution. The results for the free and total vortex 
densities are shown in Fig. |l| for A = 0.0 and 0.9. We 
do not obtain the straight line dependence of Eq. ( p^ ) 
in the lower temperature range mainly because of the 
overestimation of the free vortex number in this range. 
Nevertheless, the shape of these curves shows two slopes 
and tends to saturation in the upper temperature limit of 
these intervals, i.e. small ^'s, which is more clearly seen 
if we include more data points at higher T's. At these 
temperatures, the vortex motion is already diffusive and 
we are out of the range of ballistic motion assumed in the 
free- vortex gas theory. 

Calculations of the vortex radial pair distribution func- 
tions shows that just above Tkt the vortices are clus- 
tered with very few isolated ones (the temperatures 
T = 0.75, 0.8 for A = 0.0 and T = 0.7, 0.75 for 
A — 0.9). Increasing the temperature decreases the num- 
ber of vortices in the clusters and for the smaller ^'s in 
Figs. higher T's, the vortices are approximately ho- 
mogeneously distributed on the lattice. This beJaavior is 
similarly observed in the plane rotator modeLrnri 



B. Vortex Lifetime 

The fjcee vortex lifetime is calculated using a statistical 
methodtJ during the spin dynamics part of the simula- 
tions. The number of free vortices is counted at each 
time step dt of the time evolution and the times when 
this number decreases are saved in order to determine 
the time intervals Ati between consecutive decreases of 
the number of free vortices. Each of these intervals rep- 
resents a data point for the calculation of the free vortex 
lifetime from the following: 



lANA 



(15) 



The factor Ni/\ANi\ (Ni is the number of free vortices 
detected in the system just before the last time step and 
jAA'il is the change in the free vortex number before and 
after the last dt) accounts statistically for the possibility 
that any of the free vortices could have annihilated. The 
division by |AiVi| accounts for the case when one time 
step evolution of the system leads to a decrease of Ni by 
more than one free vortex. 

There are two requirements for the applicability of this 
method. The system has to be large enough in order to 
have a large number of free vortices for better statistics 
when using Eq. (|l5|) and the time step dt has to be much 
smaller than the characteristic time of decay Ati , which 
helps to insure that in most cases |AA^i| = 1. The time 



step dt has to be decreased when increasing T or L since 
then the number of free vortices fluctuates on a shorter 
time scale. The requirement to monitor the number of 
free vortices at each time step is the most CPU inten- 
sive task in the spin dynamics part of the simulations 
and currently makes them excessively long if we employ 
fitting to determine the precise vortex position as men- 
tioned above. We have used time steps in the interval 
1.5 X 10"" < < 1.0 X 10-3 depending on T and L and 
trying to meet the condition |AiVi| < 2. 

Each Monte Carlo IC is evolved in time up to tmax = 
100 by numerical integration of the equations of motion, 
Eq. (^), using a fourth-order Runge-Kuta scheme. Dur- 
ing the time evolution of the system, t^'s are calculated 
in order to determine an average lifetime for each IC and 
its error. The lifetimes for the different IC's and their 
errors are then additionally averaged to obtain the final 
value for the lifetime and its error for given temperature. 

The results for FM and AFM and different lattice sizes 
L are shown in Fig. |^ for A = 0.0 and in Fig. ||for A = 0.9. 
The measured free vortex lifetimes for AFM and FM, 
A = 0.0, vary between 1.03 for T = 0.75 and 0.63 for 
T = 1.1.|— Size effects are noticed for L = 32, not plot- 
ted here,E3 when T = 0.75. In this case the number of 
free vortices is very small and Eq. ( p^ ) is not reliable for 
statistical analysis. Note that the measured lifetimes for 
0.75 < T < 0.85 should be considered as an upper limit 
for the free vortex lifetime in this temperature interval 
because of the overestimation of the number of free vor- 
tices, which leads to greater values of Ni and thus of 
calculated from Eq. (|l^). The free vortex lifetimes are 
the same, within the error bars, for FM and AFM in the 
case of A = 0.0, Fig. ||, regardless of the different spin 
dynamics in both cases. The picture is quite different 
for A > Ac as we see in Fig. || for A = 0.9. The observed 
lifetimes for AFM are smaller than those for FM and this 
can be understood by the increased mobilitv of vortices 
and lower mass in AFM compared to FM.t3 

These measured lifetimes are to be compared with the 
characteristic time scale Eq. (|^) in the ideal vortex gas 
theory.B These times are listed in Table || for set of 
temperatures from the published data on the correlation 
length and the vortex rms velocity obtained from fitting 
the width and its integrated intensity of the observed 
central peak in S^^{q,uj) in simulations. The compar- 
ison of our results on the lifetime with the time scales 
in Table |l| shows that the lifetime is smaller than the 
characteristic lifetimes as much as one order of magni- 
tude for the higher temperatures listed in Table || while 
our results on the correlation length for A = 0.0 (Table |) 
are between i^i and ^2 (see Table |l|) which are obtained 
from the central peak. Even though the lifetimes we de- 
termined are shorter than the characteristic times from 
the ideal- vortex-gas theory for the temperatures listed in 
Table |l| we cannot rule out the validity of the theory to 
describe the spin dynamics because it does predict cor- 
relation lengths slightly smaller or larger, depending on 
which quantity one fits, than the correlation length we 
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determine directly and one has yet to determine the rms 
vortex velocity u independently. We cannot currently 
measure u directly from the simulations but visualisa- 
tion of the vortex positions at each time step shows that 
a free vortex almost never travels more than one lat- 
tice constant before it either becomes bound by moving 
closer to other vortices, or another vortex moves close 
to it, or vortex- antivortex pair creation occurs close to 
it. Our results on the vortex densities and clustering for 
T = Tkt + and the short vortex lifetime suggest that 
they may need be incorporated in the theory rather than 
assuming infinite vocjtex lifetime. 

Costa and Costac^ have stated that it may be neces- 
sary to consider other processes or theoretical descrip- 
tions for the cause of the CP in 5^^(q, uj), in addition to 
the ideal gas theory. In particular, they have suggested 
that vortex pair creation-annihilation events could be the 
processes responsible the CP. Since there is presently no 
theory or phenomenology that leads to any quantitative 
or even qualitative predictions, we have no way to an- 
swer this question for T > Tkt- However, for T < Tkt 
Evertz and LandauE3 have made high-precision spin dy- 
namics simulations for the FM model with A = 0. There 
they also found a CP in S^^ (q, oj) , as well as other in- 
teresting unexplained features for frequencies below the 
spin wave peak. This case is very interesting, because for 
T < Tkt the free-vortex density is approximately zero, 
whereas there can be a much larger bound vortex pair 
density. Thus we can ask whether in this case pair cre- 
ation and annihilation events may be considered to cause 
the observed CP. This question can be roughly answered 
by determining the vortex-antivortex pair lifetime Tpair ■ 
If annihilation-creation events are responsible for the ob- 
served CP, then the width of the CP should be of the 
order of Tp"„i„,. 

Our method to determine the free vortex lifetime is di- 
rectly applicable for calculation of vortex-antivortex pair 
lifetime for T < Tkt- The pair lifetime is determined 
also from Eq. (jl^) with Ni substituted by the total num- 
ber of vortices and antivortices in the system Ntot and 
|AiV,| by \ANtot\- Then Ntot/\ANtot\ will correctly be 
the number of pairs before an event of pair annihilation 
is observed, provided the time step dt is small enough 
and in a process of annihilation Ntot decreases only by 
a vortex-antivortex pair (jAA'totl = 2). We calculated 
the pair lifetime at Tkt {T = 0.7) and well below it at 
T = 0.4 with A — 0.0. The lifetimes determined are 
Tpatr = 0.487(9) and 0.39(13), respectively, with L = 64 
and 128. The average vortex pair densities for these two 
simulations are 2.4(2) x 10"^ for T = 0.4, L = 128, and 
6.9(4) X 10~3 for T = 0.7, L = 64. The pair lifetime 
for T = 0.4 may show size dependence because of the 
the very low vortex pair density for lifetime estimation 
with the linear lattice size equal to L = 128. The time 
scale of these excitations gives approximately two orders 



temperatures. This rules out the explanation of the re- 
ported central peak for T < Tkt only by the vortex pair 
creation and annihilation excitations. 



VI. CONCLUSIONS 

We carried out combined cluster Monte Carlo and spin 
dynamics simulations on classical two-dimensional easy- 
plane ferromagnets and antiferromagnets for two val- 
ues of the easy-plane strength. The correlation length 
and vortex densities were calculated in the MC simula- 
tions and their implications in the search for free vor- 
tices considered. The lifetime of free vortices and vortex- 
antivortex pairs were calculated in the spin dynamics 
simulations by averaging over the equilibrium spin con- 
figurations supplied by the Monte Carlo runs. 

The free vortex lifetime is of the same order of magni- 
tude but smaller than the characteristic time scale of the 
ideal vortex gas theory for A = 0.0 and Tkt < T < 0.85. 
This result does not rule out the validity of the theory at 
least in this interval before we have a direct way to de- 
termine the average vortex velocity. For higher tempera- 
tures and A = 0.9, the lifetime becomes smaller than the 
characteristic time by approximately one order of mag- 
nitude. The free vortex lifetime in FM is greater than in 
AFM for A = 0.9 most likely due to the greater vortex 
mass in FM compared to AFM and thus lower mobil- 
ity. For A = 0.0, the lifetimes overlap for FM and AFM 
within the error bars. 

The vortex-antivortex pair lifetimes at T = Tkt and 
T = 0.4 are approximately two orders of magnitude 
smaller than the time scale of the observed central peak 
at T < Tkt- This suggests that the pair creation and 
annihilation excitations alone can not be the reason for 
the central peak. 
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FIG. 1. The total and free vortex densities vs 1/^i for 
A = 0.0 and 0.9. The errors are smaller than the size of the 
data point symbols. 
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• 64x64, FM* 
□ 64x64, FM 
80x80, FM* 
A 100x100, FM* 
^64x64, AFM 



TABLE I. Results for the correlation length for the tem- 
peratures, lattice sizes LxL, and anisotropy parameter values 
A = 0, 0.9 considered. The correlation length is in lattice 
constant units and temperature in J/fes. 
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FIG. 2. Free vortex lifetime vs temperature for A = 0.0, 
FM, AFM, and several lattice sizes L. The FM* marked 
casep refer to previous simulations using only Metropolis et 
alEB algorithm in the MC part of the simulations. 
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TABLE IL Characteristic times = 2^i/^u, i = I, 2, 
in the ideal vortex gas theory from published data. The rms 
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